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ABSTRACT 



A simplified non-linear relay model is developed to describe ob- 
served post-stall oscillations in aircraft. The predictions of the model 
are evaluated against results obtained by numerical techniques, and shown 
to yield close agreement. 
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1. INTRODUCTION 



Considerable literature exists describing the static characteristics 
of aircraft near, or at, stall, including detailed analyses of stall 
tendencies. However, there appears to be little information on, or 
analytic treatment of, aircraft dynamics in the region around stall, and 
especially post-stall time histories in which divergent oscillations 
arise. The latter phenomenon is quite interesting, since experienced 
pilots apparently are familiar with the "rocking chair" or "porpoising" 
stall, where the aircraft, while under constant elevator control, alter- 
nately enters stall, then recovers, leading to a vertical oscillation of 
the aircraft nose about some point. 

Normally we define the airplane stall angle as that angle where the 
vehicle’s trimmed lift curve (C^ vs a) has a global maximum. It is well 
known that the lift curve then exhibits a negative slope for some range 
of a above the stall angle. Therefore, we shall also use here a^, the 
stall break angle, defined as that angle where the lift curve has maximum 
negative slope. There is a paucity of information on the nature of the lift 
curve in the post-stall or stall break region, since experimental wind 
tunnel data has usually been smoothed through this region based on engineer- 
ing intuition. In addition, an actual treatment to define the stall break 
for a physical aircraft is difficult since model data are subject to 
scaling difficulties and aircraft parameter in the non-linear region is 
limited at present. (Figure 1) 
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Typical C T vs a Curve. 



In his thesis, Frederiksen [1] considered the problem of divergent 
oscillations and limit cycles in this region. His linear analysis was 
able to show, numerically for a single aircraft, that a sufficiently 
large negative slope in in this region led to roots for the 

characteristic equation which had positive real parts (i.e. exponentially 
growing). This behavior is, of course, transitory, and disappears by 
the time the aircraft returns to a configuration where the lift slope 
is positive. However, it did provide an insight into the mechanism 
that might cause these oscillations. Furthermore, he was able to simulate 
numerically time histories in which a "rocking chair" stall did occur. 
However, due to the numerical nature of his work and his concentration 
on a single aircraft, his results are not readily extendible, or suitable 
for general analysis. 

This investigation attempts to analytically model the post-stall 
behavior in a general class of aircraft. Our main hypothesis is that, if 
the slope of the lift curve is sufficiently negative in the region beyond 
the stall break, then the aircraft will behave as if it sees a sharp, 
essentially instantaneous loss of lift. If this is the case, then the 
lift curve can be more conveniently modeled in this region by the intro- 
duction of a hypothetical "relay" acting at the stall break angle. We 
assume that, except for this nonlinearity, C^, and can be 

adequately approximated by linear functions. Those assumptions, together 
with conditions on the magnitude of various coefficients, lead to a system 
of coupled, constant coefficient, first order differential equations, with 
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a forcing term which has a jump discontinuity across a line in the phase 
plane. Analysis of this system yields a necessary and sufficient condition 
for the existence of a limit cycle, and transcendental equations predicting 
the period of this cycle. These conditions are then applied to the case 
considered by Frederiksen, and it is shown our model yields excellent 
agreement with his results. 

2. DERIVATION OF THE MODEL 

The aircraft, assumed to be in longitudinal flight, will be described 
by means of the standard coordinate system shown in figure 2. In this 
system, the equations of motion become 

{ T + qS [C sin a - C cos a ] } 

u = - g sin 0 - w 0 , (1) 



qS [C cos a + C sin a] 

w + g cos 0 + u 0 , 

m 



( 2 ) 



0 = c 

0 I C M 

y 



(3) 



and, 



. { - u sin a + w cos a} 

a = 



(4) 



(The terms in (1) - (4) are defined in the Appendix.) 

These equations are obviously non-linear, even though the coefficients 
may be linear or constant. In addition, though, in general, over wide 
ranges of variation, the coefficients are non-linear functions of the 
variables. In most applications where the range of the variables is 
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Aircraft Coordinate System 
Figure 2 




Approximation to C 

JLi 

Figure 3 
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small, they may, however, usually be adequately approximated by linear 
functions. As is obvious from figure 1, this is not the case though 

for C in the neighborhood of the stall break. 

Li 

We shall proceed under the assumption that any motion of the air- 
craft takes place in a small region about the stall break angle. Further- 
more we shall assume that and can still be approximated linearly 

in this region. We suppose the aircraft is trimmed for constant attitude 
and angle of attack, denoted by 0^ and respectively, and that the 

elevator angle and thrust, denoted by 6^ and T, are constant. The 
velocity of the aircraft in the trim configuration is denoted by v^ . 
Frederiksen's results do not show a significant variation in v , compared 
to that in a and 0 during limit cycle behavior, so we assume 



v = v 



Considering now perturbations off the trim configuration, we let 

Aa = a - a , and A0 = 0 - 0 
o o 

As noted above, we will approximate and linearly, so that 



C D " C D0 + C Da Aa » 



(5) 



and , 



C M ” C M0 C Ma Aa 2v "mS 



C I A0 



r . Aa - C M;c 6 
2v Ma M6 e 

o e 



( 6 ) 



The key ingrediant of this investigation is modeling at the stall 
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break. Reference to figure 1 indicates that, if the negative slope in 

the region past the stall break is sufficiently steep, then an aircraft 

entering this region experiences an almost instantaneous loss in lift. 

This is essentially equivalent to the action of a hypothetical "relay" 

in the lift curve function C . Thus, we shall model the lift curve as: 

Li 

C L ' °L0 + C L« 4 e + 2v~ 4 « 

e o 

+ C La Aa " C LA Sgn(a " “c } ’ (7) 

c 

where sgn(x) = |x|/x , x ^ 0 , sgn(O) = 0, and denotes the stall 

break angle. Figure 3 shows a sketch of a C curve generated by (7) 

L 

superimposed on the C curve from figure 1. (Note that, to facilitate 

L 

analysis, we have assigned algebraic signs so that, in normal cases, all 
the coefficients are positive.) 

Note now that u and w which appear in (4) can be replaced by 
their equivalent expression in terms of (1) and (2). If we do this, 
we find (since u = v cos a, and w = v sin a): 

a = - sin a - C_ + ^ cos(0 - a) + 0 (8) 

mv mv L v 



But now, using our assumption on v, and the fact that 

a = Ad , and 0 = A0 , and 0 = A0 , 



we have, 



Ad = - 



mv 



sin(a + 
o 



Aa) 



- c + - s - 

mv L v 



cos(0 



a + A0 - Aa) 
o 



A0 



(9) 
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and, 




( 10 ) 



But then, according to the assumed forms we have for C T and C„ , 

L M 

• • • • 

we see that Aa and A0 depend only on Aa, A0 , and A0 (since 

(9) can be used to eliminate the Aa term from .) Thus, for the 

remainder of our investigation, we shall be concerned only with 
equations (6), (7), (9), and (10). 

Now, under the assumption that Aa and A0 are small, we see: 

sin(a + Aa) = sin a + (cos a ) Aa , 
o o o 

and , 

cos(0 - a + A0 - Aa) = cos(0 - a ) - sin(0 - a N 

OO OO OOjAo 



+ sin(0 - a ) Aa , 
o o 



where second order and higher terms have been neglected. Using these 
relations, and the defining equation for , we can reduce (9) to 

the following: 



Ad = {- 



sm a - ^ — [C T n + C_ x 6 ] + cos(0 - a ) } 

mv o mv L0 L6 e v o o 

o o e o 



+ { 



sin(0 




cos a 



mv 



C_ } Aa 
La 
o c 



- { sin(0 - a )}A0 
v o o 



+ {1 - C 



T . } A0 + — C T . sgn(Aa -(a - a )). 

L0 mv LA 6 c o' 

o 



(ID 
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Now, using this expression to replace Aa in (7), the equation for C^, 



and substituting the resulting expression into (10), yields: 
A0 



ae - a|i [ ( - 

v L 



c „ 


- c r 6 - 


c . 


[ - — 

mv 


MO 


M6 e 2v 


Ma 




e o 




i 


mv 


<C L0 + C L6 V 


V 


cos (0 

o 



- {C + - — C . [ sin(0 - a ) cos a - C ] } Act 

Ma 2v Ma v o o mv o m La 

O O O VC 

o 



+ {%T2 sin ( 0 - a )} A0 - -r 5 - (C ' + c. [1 - -3^2 C T • ] } A0 
2v Ma o o 2v M0 Ma 2mv z L0 

o o o 



S 2 C Md C LA Sgn(Aa 



- (a - a )) 1 
c o 



( 12 ) 



Equations (11) and (12) then describe the model we shall analyze. 

In practice, many of the terms which appear in (11) and (12) are 
negligible. Calculations, using representative values, will show that 
for many, if not all cases, little is lost by assuming: 



a§c c . 


<< 1 * 


— — — c . 


2mv 2 L0 

o 




2v Ma 

o 



and also that, 



sin(0 - a ) cos a - C T 

v o o mv o mv La 

o o o c 



« 1 > 



and, 



^ sin(0 - a ) 
v o o 

o 



« 1 
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Also, we have not yet used our assumption that the aircraft was trimmed 
for flight at a Q 5 ’ This assumption says that when 



then, 



Aa = A0 = A0 = 0 , 



Aa = A0 = 0. 



Thus, imposition of the trim conditions reduces (11) and (12) respectively 
to : 

0 = { — sin a - [C A + C T x 6 ] + cos(0 - a ) } 

mv o mv L0 Lo e v o o 

o o e o 



- C T A sgn(a - a ) 
mv LA ° c o 
o 



and, 



0 = “ C ^n " C ma 6 
MO Mo e 

e 



(13) 



(14) 



If we now incorporate (13) , (14) , and our assumptions about small 
terms into (11) and (12), we arrive at 

Ad = A0 + C T . {sgn(Aa - (a - a )) + sgn(a - a )} , 
mvLA & c o ° c o 

o 



(15) 



and. 



A6 = 



- C Aa 
I • Ma 

y 



-2 

qSc 

2v I 
o y 



(C • + n . } A0 
M0 Ma 



- 0,. C {sgn(Aa - (a - a ) + sgn(a - a )} . (16) 

2ml v^ Ma LA 6 c o & c o 

y o 



10 



We now define the quantities 3> a), a, and £ by: 



I Ma 

y 



(6 2 + w 2 ) > 0 



* 



qSc 2 

2v I 
o y 



{C. 



Me 



+ C M .} = 26 > 0 
Ma 



( 17 ) 



mv 

o 



LA 



= a > 0 , 



and. 



2v I Ma 
o y 



C > 0 . 



Since it becomes important in our later discussion, note here that 
the quantity: 



-2 



(2B - » - fri c h9 ‘ 0 • 



o y 

Thus our model for the flight dynamics of an aircraft, trimmed for 
longitudinal flight at an angle of attack near the stall break angle, 
becomes : 



Aa = A0 + a(sgn(Aa-(a - a )) + sgn(a -a ) } . 

& c o G c o 

• • 2 2 

A0 = - (6 +03 )Aa - 2 3 A0 - £a{sgn(Aa-(a - a ) ) + sgn(a -a )}. 

co co 



(Note that (18) and (19) describe a system which is second order in Aa 
and A0 . Since the original system (1) - (4) was equivalent to a fourth 
order system in a, 0, 0, and v, our model will obviously not describe the 



(18) 

(19) 
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full range of dynamics. However, we will show that it does appear to 
predict the known limit cycle behavior. In essence, what we are doing 
in equations (18) and (19) is basing our analysis on the influence of 
the stall break on the short period mode of the aircraft, and neglecting 
the long period mode. Furthermore, the mathematical advantages of dealing 
with a second order system are enormous.) 

3. FLIGHT WITH TRIM AT THE STALL BREAK ANGLE 

We begin our investigation of the model described by (18) and (19) 

with the case a = a , that is the aircraft is trimmed for flight at 
co 

the stall break. Then since a - a =0, the second sgn function in 

co & 

both (18) and (19) disappears. If we let = A0 and = Aa , our 
system reduces to: 

2 2 

Xi = - 26 x 1 - (6 + co ) x 2 - Ca sgn 

x^ = x^ + a sgn x 2 (20) 

A. The Phase Plane . Consider the analytic behavior of this system 
in the (x^ , x^) phase plane. The line x^ = 0 divides the plane into 
two regions, where separately the equations are linear, constant 
coefficient and non-homogeneous . For example, in x^ > 0 , the system is: 

2 2 

x^ = - 20 x 1 - (6 +o) )x 2 - £a , 

x 2 = x 1 + a . (21) 



12 



A key observation now is that (20) is symmetric, i.e. 



{ (x i (t), X 2 (t)) I T 1 < t < T 2 } 

represents an arc of a trajectory if and only if 

{ (-x 1 (t) > -x 2 (t)) I T 1 £ t < t 2 } 

also does. Thus we need direct our attention solely to the region 
> 0, and draw any needed information about < 0 from symmetry 

arguments. We shall use the usual convention in referring to the line 
x^ = 0 as the switching line. 

We start by considering the critical points of the system, i.e., 
the points where x^ = x^ = 0 . These are given by: 

X 1 = X 2 = 0 » 

(23 - C) 

X 1 = a ’ x 2 = ' W +” V • ’ 

and, 

(26 - O 

X 1 = _a ’ X 2 = (6 2 + A a 

(These critical points actually occur, since, as we have noted, 

(26 - C) > 0 ). The presence of the sgn term in (20) implies that the 
critical point at (0,0) will be unstable. However, it is easily shown 
that the homogeneous part of (21) has the characteristic polynomial 

A 2 + 26A + (8 2 + a 2 ) , 

whose roots must have negative real parts according to our definition 
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of 3 and a) in (17). Hence, the two critical points not at the origin 
are locally asymptotically stable. (That is, small perturbations off 
this point decay back into it.) We shall show shortly these are not 
globally asymptotically stable. 

B. Conditions for a Limit Cycle . We now consider the conditions 
under which (20) can have a limit cycle (i.e. periodic) solution. 

Obviously from the stability discussion above, any limit cycle trajectory 
must pass through the switching line. Furthermore, it must by symmetric 
with respect to the origin. Thus, if we consider the arc of the limit 
cycle in the region > 0, and assume this arc originates on the 

switching line at the point (x*, 0), then, after a half-period, it 
must terminate at the point (-x*, 0) . Thus, if we denote the half-period 
by T , we have a necessary condition for a limit cycle of period 2T is: 

x x (0) = - x x (T) , 

( 

x 2 (0) = x 2 (T) = 0 . 

(This situation is represented in figure 4.) 

It is easily shown that in the region x 2 > 0 the solution to 
(20) is given by: 





(23) 



and 



* 2 (t) = c j_ e 



-6t 



cos tut + c 




(24) 
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The Limit Cycle Necessary Condition 
Figure 4 




The > 0 Region 
Figure 5 
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(These simply represent arcs of counterclockwise spirals centered at the 



point (-a, + ^2 ) a ) > which is the critical point in the upper half 

plane. Thus every curve originating on the switching line must follow 
a trajectory that carries it above this point. (A typical such curve is 
labeled (1) in figure 5.) Furthermore, since the critical point (if it 
exists) is asymptotically stable, there must be some curves entering the 
region which never exit it, but decay down to the critical point. (One 
of these is shown as ( 2 ) in figure 5.) Finally, since any trajectory 
intersecting the line x^ = -a has a horizontal slope at the point of 
intersection (see eqn (21)), then the boundary between the decaying curves 
and those that exit to x^ < 0 must just be that curve which terminates 
at (-a,0) . (This curve is labeled (3) in figure 5.) 

But also note that the time interval represented by a trajectory in 
the region x^ > 0 is a function of the counterclockwise angle (as 

(2R — r ^ 

seen from (-a, Z\ a ) ) between the initial point (x 1 (0) , 0), 

Cp + co ) -L 

and the point at which the trajectory first touches the switching line 

again. Thus, any curve which represents a half-period of a limit cycle 

must spend less time in x n > 0 than T , where T is the time 
r 2 o o 

represented by the curve (3) in figure 5. Thus T satisfies: 

T < T , 
o 

where , 

x-(T ) = -a 
1 o 

x (0) = x (T ) = 0 . (25) 

z z o 
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Conditions (22) and (25) can now be seen to provide both the necessary 
and sufficient conditions for existence of a limit cycle of period 2T . 

The condition x 2 ^) = x^(T) = 0 , applied to (24), immediately 
leads to: 



and , 




(26 - c) 

7F"+ a ’ 

cos wT) = c 2 sin toT . 



(26) 



(27) 



Notice, according to these, that T ^ — . 

CJ 

The second condition, on x^(t), when applied to (23), implies: 



x^(0) “-Bc^+wc^-a"- x^ (T ) 



— RT 

= c^e (6 cos rnT + m sin mT) 

_ 

- C 2 ^ (6 cos coT — to sin toT) + a 



or. 



— ftT — RT 

c^(B + Be cos (joT + me sin mT) 



— RT — RT 

- <^ 2 ^ + 006 cos wT ” sin wT ) = “ 2a 



(28) 



However, we have already commented that 

TT 

T ^ — (because (26 - £) ^ 0) . Then sin u)T ^ 0, and hence we can multiply 
m 

RT 

both sides of this last equation by the term e p sin mT, and simplify 
{using (27) to replace the term c^ sin mT) to yield: 

c 1 (26e 0T sin wT + oj - ue 20T ) = - 2ae BT sin wT , 



or, 



sinh 6T 



(2a + Bcp 



OJC 



1 



sin uT , 
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which, using (26), reduces to: 

= inh st - - ( » 2 o b2) sin » t • < 29 > 

Equation (29) expresses the necessary condition. In order to guarantee 
sufficiency, i.e., in order that a solution of (29) must generate a 
limit cycle, we must add the restriction, 



T < T , 
o 

where T was defined above, 
o 

Applying the definition of T q to the solutions given by (23) 
and (24) yields: 



and , 



„ = <23 ~ O 

1 (8 Z + o» z ) 



c, (e^° - cos ) = c- sin o>T , 



- 0,6 (8 cos wT + sin aiT ) 

1 o o 



_ gT 

+ c 0 e 0 (to cos ojT 
2 o 



- B sin ojT ) = 0 
o 



Note that the first two equations are identical to (26) and (27) . This 

follows since the conditions on X£(t) implied by both (22) and (25) 

are identical. Again we can assume sin ajT^ / 0, multiply both sides of 

RT 

the last equation by e P ° sin simplify by roplacing the term 

c^ sin ayr^ by its value in terms of the middle equation, and divide by 
c^, to yield; 

_ g T 

(o cos ooT - 8 sin wT = co e ° 
o o 



(30) 
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In view of our earlier comments then, we must pick the smallest 
positive solution of (30). 

Mathematically, then we can summarize our results in the following: 
Theorem : The system (20) will have a limit cycle of period 2T if and 

only if T satisfies: 

2 2 

sinh BT = - B sin o»T, (29) 

o) (26 - c) 

and , 

T < T , 
o 

where : 



T is the smallest positive solution of 
o 

_ grji 

o) cos o)T - 6 sin o)T = w e o . 
o o 



(30) 



C. Example - The Case Considered by Frederiksen [1] . To 
determine validity of the predictions of this model, consider the case 
used by Frederiksen [1]. Figure 6 shows the table of values appropriate 
to this case, which describes the F-94, a straight wing, single engine 
jet. The value of the coefficients appropriate to the linearized 
differential equations were computed and are shown in figure 7. Note that 
our assumptions on small size of certain coefficients seems valid. Using 
the values of 6, a) and £ so computed, equations (29) and (30) were 
graphically and numerically solved, and are graphically portrayed in 
figures 8 and 9, respectively. There are obviously two solutions of 

7T 2 IT 

(29) in the region — < T < — (which is the region of interest since 
(26— ^ ) >0.) These are 



T = 1.87 and T = 3.56 
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p = 2.38 x 10 
m = 384 
S = 239.0 



|1 

Parameters for F-94 



c = 6.40 
I = 2.65 x 10 4 

y 



C L0 mlA *° 



C T = 0.30 
La 

c 



C T . = 0.05 
LA 



°M0 



= 0.15 



C„ = 1-54 
Ma 



% - 8 - 16 



0 - 0.099 rad 

o 



a =0.41 rad 
o 



a =0.41 rad 
c 



jl_ Values used by Frederiksen [1] 



C L 6 ' °' 43 
e 




Sla 



4.27 



C M6 =°- 88 
e 



v = 169 ft/sec 
o 

T = 0.0 



Figure 6 
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Based on the parameters in figure 6, 

qS° o 'T 

C T • = 2.4x10 , -S- r = 0 .O 8 I 

2tnv 2 L0 2v M “ 

o o 



sin (6 - a ) = - .058 

v o o 

o 



sm (0 - a ) 

v o o mv 

o o 



cos a - — C T 

o mv La 
o c 



= - .096 



(jSc ^ 2 

— C Ma " (f? + “ > = 3 ‘ 00 

y 

—2 
qSc 

2v I 

o y 



{C M0 + C Ma } = 2B = °- 46 ’ 30 6 



= 0.23, u = 1.72 



-3 

C T = a = 6.3 x 10 

mv LA 
o 



Sc 



2v I Mi 
o y 



c,, = C = 0.16 



(2g-c) = 0.30 > 0. 



Coefficients for the Differential Equations 



Figure 7 
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Graphical Solution of (29) 
Figure 8 




Graphical Solution of (30) 
Figure 9 
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But T , defined as the first positive intersection of the curves 
in figure 9, is given by 

T =3.10 
o 

Thus, only the lower value of T satisfies the requirement 

77 

— < T < T and there will be exactly one limit cycle, with period 
03 o 

2T = 3.74 sec 

Computing c^ and c^ according to (26) and (27) , and using the 

values to describe X£(t) = Aa(t) as given by (24), leads to an 

amplitude for this limit cycle of 

Aa = .018 rad = 1° 

max 

These values are in very good agreement with those obtained by 

Frederiksen, which were 2T = 3.77 sec 

Aa = .020 rad 

max 

D. Additional Comment . We shall make one other fairly important 
observation on the mathematical nature of (29) and (30) . Note that (29) 
is equivalent to 

sinh gT _ _ 03 2 + BC ~ B 2 
sin coT 0)(26-C) 



Now since (2B~C) > 0, the region of interest for this equation is 

TT 2tT 

— < t < — . Letting 

03 03 & 



G(t) 



sinh Bt 
sin 03 1 



(31) 



(32) 
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we have 



and, 



Furthermore 



G(t) < 0, in ~ < T < — 

03 03 



* 



< + > = " - 



co 



G'(t) 



B cosh Bt sin 03t - 03 sinh Bt cos 03t 
. 2 

sm 03 1 



Note the presence of the square in the denominator implies that the 
algebraic sign of G f (t) is the same as the sign of the numerator. 
This numerator, 

B cosh Bt sin 03t - 03 sinBt cos oit 

has the value 

• -i B^r ^ r\ _ _ tt 

03 sinh — > 0 at t = — , 

03 03 

and the value 

. i 2 Btt „ 2tt 

-03 sinh < 0 at t = — 

03 03 

The derivative of this numerator is 

2 2 

(3 + 03 ) sinh Bt sin 03t. 

Since this derivative is negative on ^ < t < ^ , then it clearly 

7T 

follows that the numerator of G t (T) starts positive at — , remains 
positive up to some point T , — < T < , then turns negative, 

r r r m 03 m 03 

staying negative for the remainder of the interval. Based on earlier 
comments, G T (t) behaves the same. Thus G(t) is monotonically 
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7T 

increasing on — < T < T m , and mono tonically decreasing on 

2tt tt 2 tt 

T < T < — . Hence G(t) has a unique maximum in — < t < — 

m a) n a) o) 

[Figure 10 shows the general form of G ( t ) . ] 



But now let’s investigate the relative positions of T and T 

o m 

At T , G’(T ) has the numerator 
o o 



B cosh BT sin o)T - o) sinh BT cos o)T 
o o o o 

1 BT 1 _ bt 

= — e P o [B sin o)T - o) cos o)T ] + — e ° [B sin coT + o) cos o)T ] . 

2 o o 2 o o 

But according to (30) , 

_ 

B sin 03 T = - u)e o + o) cos o)T 
o o 



Thus, the numerator of G’(T ) is 



1 BT 0 r -BT . 1 -BT - T 0 , 0 m 1 

= ^ e ° [-we o] + — e o [-o)e + 2u)cos 



= ~ e BT o [cos u)T - cosh o)T ] < 0 
/ o o 



Then 



G’(T ) < 0 , i.e,, G(T) is decreasing at T , and so by 



our earlier comment, T < T . 

m o 



But observe that the right hand side of (31) is independent of T . 
Thus, since (29) and (31) are equivalent, a necessary and sufficient 
condition for (29) to have a solution, based on our analysis of G(T) , is 



max 

TT . 2TT 

- < t < 

0 ) 03 




2 2 
> + Be - 8 

oo(2B-C) 



(33) 
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Figure 10 




Figure 11 
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7T 4“ 

Finally, since G(— ) = - 00 and T < T , then if (33) is satisfied, 

a) mo * 

then there must be a T < T < T such that 

m o 



2 2 

sinh 3T _ _ co +3^-3 
sin o)T o)(23-£) 



Then this T satisfies (29) and (30), and so generates a limit cycle. 
Thus, (33) express a necessary and sufficient condition for existence of 
a limit cycle in this model. 

This is expressed by the following: 

Theorem : The system (20), where (2g-£) > 0 , will have a limit cycle 

solution if and only if 



max 

7T . 2„ 

— < t < 

0) (A) 



[ sinh gt 
sin wt J 



> 



2 2 
a) + - B 

w(20-O 



Observe that, with the substitution 



A = — > 0 , and rj = ~ > 0 , 

(JL) 0) 



the necessary and sufficient condition becomes 

2 



max 

7r<t<2.TT 



[ sinh At 1 > 1_ 

sin t 



+ T)\ - A 

(2A - n) 



(34) 



(the left hand side of which is clearly a function of A only.) Note 
we already have commented that (2B - ?) > 0, hence 

(2A - n) > 0, or 0 < ri < 2A 
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Now, observe that 



d_ 1 + lU - X 2 1 _ x 2 + 1 

dn [ (2X - n) J (2X - n) 



< o 



and therefore the right-hand side of (34) has its maximum value at 
C = 0. Thus, clearly, a sufficient condition for (34) to be satisfied, 
and hence for a limit cycle to exist, is, 

max 
7T< t < 2 IT 

In figure (11) we plot the function 
max 

7T<t<2Tr 

for 0 < A < 1, and observe that (35) cannot be satisfied if 



sinh At j 1 - A^ 

sint J 2A 



sinh Xt 
sint 



X 2 - 1 

2A 



(35) 



X > .280 . 



(Observe (35) cannot be satisfied for A > 1.) Thus, a sufficient 
condition for existence of a limit cycle is: 

X < .280 (36) 



Now, if A< .280 , and 8 = Ato , 
we have: 

9 ? 2 2 2 2 
to < to + B = (1 + A ) to < 1.078 to , 

[~2 2~ 

and so we can approximate, to = v to 4- 3 

But then, according to the definitions (17) , X = ~ is proportional 
to 

C . + C ’ 

Ma M9 

> 

C Mct 
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and we see our sufficient condition (35) translates into a statement 
that 

C . + C * 

Ma M0 

is somehow small compared to C w 

r Ma 

Since is analogous to a spring constant, and C^* and are 

damping terms, we can view (35) and (36) as saying: 

M A sufficient condition for existence of a limit cycle behavior 

in the presence of a sharp negative lift slope, is that the damping 

terms in C , are small relative to the "restoring" term, C 

M Ma 

Clearly this interpretation makes physical "sense". 

4. SUMMARY AND CONCLUSIONS 

In this paper, we have developed and analyzed a simple two-dimensional 
model for post-stall time histories in the presence of a sharp lift loss 
at stall. This model yields a simple necessary and sufficient condition 
for the occurrence of limit cycle, i.e., "rocking chair" stall. Further- 
more, the model predicts values that are in very close agreement with 
those obtained from numerical solution of the full, four-dimensional, non- 
linear case. Lastly, the model yields physically satisfying information 
about the important parameters for occurrence of this limit cycle. 
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APPENDIX. Definition of Variables and Coefficients 



c = Mean Aerodynamic Chord (MAC) 

C D = Drag Coefficient 

C T = Lift Coefficient 

C^ = Pitching Moment 

g = Gravitation Acceleration 

I = Moment of Inertia about c.g. 

y 

m = Mass 

x 2 

q = % P v 

S = Wing Reference Area 

T = Thrust 

u = Velocity Component along Thrust Line 

v = Velocity 

w = Velocity Compenent normal to Thrust Line 

a = Angle of Attack 

8 = Pitch Angle 

p = Density of Air 
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